Electron gas under the non-homogeneous magnetic field of cylindrical symmetry 
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We study the states of SD-electron gas in non-homogeneous magnetic field. It is supposed that the 
step of magnetic field, at which field changes its sign, Hes on the cyhndrical surface. The eigen value 
problem is solved for are different parameters of the system. The equilibrium statistics in the limit 
of large quantum lifetime are considered as well. It is predicted that the system characterized by 
quasi-zero-dimensional spectra with large density of states; snakes orbits lead to the effective charge 
transfer; an appearance of the magnetic field step strongly modify the coordinate dependence of 
carriers concentration. 

PACS numbers: 



I. INTRODUCTION 

The fabrication technologies for making non- 
homogeneous magnetic field in 2DEG are well developed. 
Micromagnets, superconducting micro samples provid- 
ing Meissner effect, or magnetic impurities are used 
to obtain non-homogeneous magnetic field with given 
strength and geometry [TJ |5] . 

Different phenomena related to 2D-electron transport 
under non-homogeneous perpendicular magnetic field 
have been studied intensively within recent years both 
experimentally [2HB] and theoretically [3 [3 [TOI? 41^. 
If the characteristic geometry scales of non-homogeneous 
magnetic field Ich satisfy the condition Ich > le (le is 
the electron mean free path), it leads to the drift diffu- 
sion phenomena. If the contrary condition l^h < h is 
satisfied, but the temperature T or characteristic relax- 
ation energy ^ h^{llfh)~^ ~ '''qi^^ of the system is large 
enough than classical ballistic transport takes place 
(Tqi is the quantum life time, fh is the effective mass of 
electron). The quantum phenomena takes place if the 
average of absolute value of cyclotron frequency satisfy 

\^\>Tn-\T^^\ 

In this paper we propose to examine properties of SD- 
electron gas in non-homogeneous magnetic field in quan- 
tum and semiclassical ballistic regimes uJc 

The ways to obtain 3D-electron gas under strongly 
non-homogeneous magnetic field are not widely discussed 
to the best of our knowledge. The appearing of the high 
density currents which are arranged in are special way 
may cause strongly non-homogeneous magnetic field of 
requisite geometry. 

The superconducting inclusions in bulk semiconduc- 
tors have been examined for several decades. More- 
over, present-day nanotechnology permit to obtain hight- 
quality superconducting nanowires of given geometric pa- 
rameters, using different superconducting materials and 
substrates |l3lll4j. But most of the studies that had been 
carried out, aimed to study superconducting phase tran- 
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sition, and tunneling effects (proximity) of the system 
under different parameters and conditions. There are 
number of conditions which could be responsible for the 
phase transition in the real system with superconducting 
sample (sample geometry, temperature, magnetic field, 
irradiation etc.) |15l[5D] . One of the most significant pa- 
rameters for superconducting nanowires is their diame- 
ter. The decreasing of the diameter of nanowire to the 
order of superconducting coherence length or lower may 
dramatically change or completely vanish the supercon- 
ducting properties, more over quantum confined effects 
are also possible [2Tti25] . 

Specially arranged currents integrated into a bulk 
semiconductor may be the source of strong non- 
homogeneous magnetic field with given geometry. The 
examples of fiat and cylindrical geometry of the magnetic 
field step are shown in Fig. 1(a) and (b) correspond- 
ingly. One can use classical limit of the model to make it 
clear. In classical limit far from rupture of magnetic field 
charged particle is in homogeneous magnetic field and 
has spiral-like trajectory with parallel to the magnetic 
field axis. On the other hand, if the particle closer to the 
currents than its cyclotron radius, the trajectory, which 
is schematically described in Fig 1(d), become possible 
(so-called snake orbit) [8]. Obviously, the flat geometry 
case (a) for the most part is similar to the studied in 
[SI H] case of perpendicular to the 2D-electron gas non- 
homogeneous magnetic field. Therefor the second case 
[Fig. 1 (b-d)] is more interesting 

The currents /i, /2, are arranged on the cylin- 

drical surface and of the same directions and values. The 
current /q = — (X)fc=i -^fc)/^ (opposite direction) is fixed 
to the axis of the cylindrical surface. 

We make the next simplifications. The semiconduc- 
tor electrons has parabolic isotropic dispersion law. Any 
reasons to proximity effect appearing (contact peculiar- 
ities etc.) are absent. The Fermi energy of the electron 
gas is much lower then threshold of superconductor. The 
superconducting inclusion with /q (see Fig. 1(d)) may be 
treated as impermeable barrier. The distance between 
the currents on cylindrical surface /i, I2,...,In and the 
diameters of the corresponding inclusions is lower or in 
order of particle wave length [Ich ^ Ae). We also suppose 
that the system may be treated as infinite one and has 
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FIG. 1: Schematic description of the systems with magnetic 
field step; a - plane arranged parallel currents, and possible 
trajectory of the charged (negatively) particle; b - currents 
are arranged on the cylindrical surface and one fixed to the 
its axis; c - coordinate dependence of the magnetic field and 
vector potential (arbitrary units) ; d - the section of the de- 
scribed in (b) geometry and possible trajectory of the motion 
of the particle. 



continues translational symmetry along the the cylinder 
axis. The spin effects are not considered as well. 

Also we are not considering any mechanisms of relax- 
ation which cause electron spectra transformation and 
broadening of the energy levels (limit of large quantum 
lifetime) . 



II. SCHRODINGER EQUATION 



form = exp(-iPz - im0)$"''='/'=-(rfcc). Before 

substituting this form into eq. (1) we turn to the next 
dimensionless parameters, wave number = /k^, ra- 
dial coordinate p = rk^ [pc = rckc, Psc = fsckc see Fig. 
1(d)]. Substituting definite above form, one obtains next 
homogeneous differential equations for 



1 d d vn? 
p dp^ dp p^ 



In 



-In 



P 



0, 



(2) 



where k} — 2mEm,k^ I {iT-k, 



- The bound- 

ary condition (i) - ^pi!^'^' = is the sequence of 
impermeable barrier approach (see above), other one 
is in standard requirements. The standard require- 
ment for such solutions are next (ii) - the existence 
of the integral \^p'"^''^^\^pdp, in particular expects 

for p\^p'™''^ Ip-i-oo ^ 0; (iii) ^ normalization condition 
/o°° l^p'"^"'^^ l^pdp = 1, always may be satisfied if (ii) is 
satisfied, but exclude one constant from the general solu- 
tion of eq.(2). (iv) - the continuity of the first derivative 
+ o)/dp - d^'^'^^p - o)/dp]o^o = (or ex- 
istence of the second derivative) for p € (psci oo); The 
condition of the orthogonality of different eigenstates is 
also satisfied. 

It is convenient to make substitution x ~ ln{r/rc) that 
results to the usual form of one-dimensional Schrodinger 
equation 



d^ 
dx^ 



0, 



(3) 



Now we turn to the quantum mechanical description of 
the problem. The Schrodinger equation for one-electron 
wave function in approach of effective mass has the next 
form [? 1: 



2m 



E„ 



(1) 



here fh is the effective mass of electron, e = — |e| is 
the electron charge, c is the light velocity, Ar is the 
vector-potential. The wave function that describes 
the eigenstate with with eigenvalues of the quantum 
numbers k^, m, E„i^k^. One has choose the vector- 
potential to not break continuous translational symme- 
try along cylindrical axis, for example Ar = (0, 0, A^^) 
with nonzero z-component = {pcc/e)\ ln(r/rc)|, where 
Pc = fikc — 2e/o/c'^ [here and below we use cylindrical 
coordinates r = (r, (/), z), all vector components has 
corresponding upper indexes]. Notice, that we are not 
considering field inside superconductors, and so gauge 
invariance has trivial form. The momentum operators 
— eA^ /c and are not commuting with hamiltonian, 
but and rp* = —ihd/d(j) are commuting. The lat- 
ter implies that we can search the solutions in the next 



were if^'™'''^ = exp(2a;)/9^ [k? -t- 2k^\x\ — x^) — (in- 
troducing one-dimensional potential energy U(x), one 

obtains if^.™."^^ = ^[Em,k- - U{x)]) Below we ana- 
lyze this equation using quasiclassical approach. First, 
we suppose condition X2 — Xi ^ 1, were xi^2 are turn- 
ing points: K!^i^'^^ = 0. The solutions inside the region 
xi <^ X X2 tend to the functions [57] 



(4) 



where C+ and C_ are constants. The condition (i) trans- 
forms to = [xsc = ln(rsc/'^c)]- The solutions 
for X2 ^ X approaches the functions [conditions (ii)]: 



(][)« 



i3exp{-/;/Cy"--^dx'} 

V l^x 



(5) 



where /C^'™'""^ ^ m? - exp(2x)/9^ (k^ + 2K^|a;| - x"^). 
The condition (ii) for approach (4)-(5) means that if x 
{x < X2) is increasing and passing by point X2 in complex 
space {x ^ x' = X2 + 5xe^'^^ , 5x =\x — X2\> o) then (4) 
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FIG. 2: Spectra of the system versus for different Jq: (a) 
Jo correspond to psc = 4:, pc — 10 and (b) Jo correspond to 
Psc = 2, Pc — 5. 



turns into (5) [2T. There are three kinds of quasiclas- 
sical solutions which satisfy (i)-(iv) and correspond to a 
different choice of C±, k, in (4). 

The first one (k^ ^ 0, <|^.'".'^^ - antibonding states): 



sm 



A- 



(6) 



with = 7r(n2 + f ) and £^ K^r^'^' dx' = 

7r(ni + 1) if Xsc > xi, or j"^ K^;"^''' dx' = 7r(ni + |) if 

Xsc ^ X\. 

The second one (k^ ^ 0, (j)^'™."^ - bonding states): 



cos ( Jo" K:"'^'' dx' 



V "a; 



(7) 



with if^";™^'^^^^' = 7r(n2 + i) and £^ i^«;'"^''^dT' = 
7r(ni + 1) if Xsc > xi, or j"^ iiT^;™'" dx' = 7r(ni + i) if 
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FIG. 3; The radial wave functions versus p for pc = 10, 
Psc ~ 4 and different numbers ni, n2 or n [see eqs. (6), 
(7) or (8)] and m: (a) antibonding states (nl, n2, m) = 
(0, 2, 0), (0, 1, 1), (0, 3, 0), (0, 3, 4); (b) bonding states 
(nl, n2, m) = (0, 2, 0), (0, 2, 4), (0, 1, 3), (0, 1, 4); states 
with = [n, m) = (1, 0), (3, 0), (0, 2), (0, 3); 



Xsc < x\. Here, in eq. (7) we suppose > +o (see 
Appendix B and [? ]). 
The third one (k^ = 0): 



sm 



(/.I KT^^^-'dx^ 



(8) 



with J^^ iiT^;'"'" °c?x' = 7r(rt + |) if Xsc < xi and 

SxL^x'"^'" = ■^(" + 1) if Xsc < xi {n, ni, 7i2 

are integer numbers). 

Using (5)- (7) we construct the solutions of the problem 
(2), (i)-(iv) in the first approximation. Next approxima- 
tions performed numerically. In Fig. 2(a),(b) one can 
see the spectra of the system versus quantum number 
k^. The numerical solutions which correspond to the 
first approximations (5)-(7) are shown in Fig. 3(a)-(c) 
correspondingly. 

All states with m 7^ are twice degenerated E^m^ j^z = 
The lager the lower difference i?|m|+i,fc^ ^ 



E_ 



|m| 



> 0. The bonding states and the states with 



= may be strongly modified or suppressed because 

Pa 



(j,K,m,K _^ Q g^^j they interact with the currents on cylin- 



drical surface r = r^. On the other hand the states with 
= plays an essential role in intersubband photoexci- 
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tation processes if the Fermi energy of electron gas and 
photoexcitation energy Ep, huj < p^/2m. 



III. EQUILIBRIUM DENSITY MATRIX 

Now we can turn to the consideration of the statistics 
of the system under equiUbrium conditions. We neglect 
any effect of the scattering mechanisms on eigenstates 
and eigenvalues of the system, supposing that feeble in- 
teraction between electrons and thermostat, as well as 
electron-electron interaction impose equilibrium state of 
the system that may be characterized by diagonal density 
matrix: 



\E,k^,m){E',k'\m'\ 
cxp[{E^,k^ ~ fi)/T] + 1 



(9) 



here Sa^t is the Kronecker symbol, fi and T is the chemi- 
cal potential and the temperature of the system. The the 
coordinate-dependent concentration is calculated with 
the usual formula 



N(r) = 2Spg, 



(10) 
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where SpA denotes the sum of diagonal matrix element, 
multiplier 2 takes into account spin. In Fig. 4 one can see 
the concentration versus dimensionless radial coordinate 
p calculated using eq. (10) for different Iq (variation of 
Iq means proportional variation of pc and psc) P and T. 
Comparing Fig. 4 (c) and (d) show the variation of the 
concentration with magnetic field step appearing. In the 
case lo ^ 0, J2k=i h = and > 10 if T < pl/2m 
the majority of electrons localized inside the cylinder 
r < Tc- Classically they are not effected by the mag- 
netic field modification after switching on the currents 
in the cylindrical surface -^fc 7^ 0' t>ut their coor- 

dinate distribution are modified essentially. The sketchy 
description of the effect one can see in Fig. 5. 

Using the current operator = e{pz—pc \ ln(r/rc)|)/m, 
we obtain 



= 2Sp{fQ). 



(11) 



One can see the current density versus p for different T 
that is calculated using eq.(ll) in Fig. 6. The current is 
increasing with the increasing of p, T, and Tc as well. 
The the reason of the increasing are the same versus p, 
T, and Tc - the higher one of this parameters the more 
states participate in charge transfer. The increasing of 
Tc leads to increasing of the density of states. And the 
increasing of T and p leads to enabling of higher states. 



IV. CONCLUDING REMARKS 

To summarize, we have considered electronic gas in 
discontinues magnetic field under cylindrical geometry 



FIG. 4: Concentration versus p: (a) T = 0, /i = 0.8pl/m for 
Psc ~ 4 and pc = 10, the lines 1-3 take into account antibond- 
ing, bonding and the states with — 0, the line 4 includes all 
of them; (b) the same as in panel (a) for p^c = 2 and Pc = 5; 
(c) r = and different chemical potential value p = 0.8 x, 
1.6x, 2.4x, 3.2x pl/2m marked 1-4 correspondingly; (d) the 
same as in panel (c) for the case Jo 7^ 0, X/fc!=i Ik ~ (cur- 
rents on cylinder are switched off) [compere (c) and (d)]. 



condition. The quasi-zero-dimensional spectrum was ob- 
tained for semiconductor electronic gas under consid- 
ered condition. Nevertheless average value of the cur- 
rent along the preferential direction is predicted to be 
non-zero. 

Now we discuss the numerical parameters that where 
used. For example, one can suppose that system based 
on bulk GaAs and lead nanowires. Under the helium 
temperatures characteristic thermal momentum of GaAs 
electronic gas pt — VSmT ~ 2.7 x 10^^^ g ■ cm/ s. Sup- 
posing — 0.22 pm we obtain the Ginzburg - Landau 
parameter for such lead nanowire under T ~ 4.2K es- 
timated to be 0.31 < l/-\/2, so we deal with first type 
superconductor (see also [25^ ) . The critical field of lead 
under helium temperature He — 530 Gs. Using the Sils- 
bees rule one can obtain pc ^ lApT and VscPc/^ — 4. 

Finally, we estimate the concentration and current val- 
ues. Using estimated here above parameters one can 
make conclusion that fcj? ~ 5 x 10^^ cm"'^. 

We performed consideration suppose very hight cur- 
rent densities /q, Ii,...In (i-e. low psc and pc) to em- 
phasize quantum peculiarities. The realization of such 
systems is questionable, but the majority of discussed 
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FIG. 5: The sketchy description of the magnetic field step 
effect on the electrons coordinate distribution. The effective 
one-dimensional potential energy 1]q{x) — h^{k^'^ + 2k^kcX — 
k^x'^)/2m correspond to the case X/fc-i = 0, fc^ < [see 
eq. (3)]. After the currents are switched on X^^j^ Ik 7^ the 
effective potential energy is changed to the U{x) = Ti^{k^^ + 
2k^kc\x\ — k1x'^) /2m (k^ is the same). The electrons tunneling 
to the right hand side (outside the cylinder). 

peculiarities remain unique even for the case of low cur- 
rent density [i.e. 1 ^ psc ^ Pc, see appendix A and the 
discussion of Fig. 4(c), (d), Fig. 5]. 
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Appendix A: Flat geometry as the limit case 

Here we consider eq. (1) in the limit case Vc ^ 00. 
We locate the Cartesian coordinates in the point (f> = 0, 
r — Vc, z = and make asymptotical transition from 
the cylindrical to the orthogonal coordinates x = r — Tc, 
y = rcf/), z = z {4> posses low values ~ r~^, for finite y). 
Making expansion of the vector-potential up to second 
order on x/rc and substituting wave function ^(r) — 
^K,,K ii^-j gxp(«fcyy -I- ik^z) one could obtain 



dx2 



+ 2k^k, 



2r2 



2 ~ 2 



(N-xo)^ 



(x) = (Al) 



here fc^ = 2'mE/fi^ — fc^^, xq = Vck^/kc = P7.l{fhujc), 
and ujc — Pc/ifcfn)- This equation still do not corre- 
spond to flat geometry case Fig.l (a). One has to use 
the next additional restrictions: (i) ^ kc — cjcfcfn/h, 
(ii) rcexp{—k^ /kc — k/kc) > Tsc to obtain flat geome- 
try case. The restriction (i) we need to neglect the term 
2fc^fcc|x|x/(2r^). This term means that system remem- 
bers cylindrical symmetry even for large Tc if quantum 
number fc^ is large enough. The restriction (ii) we need 
to omit boundary condition for r = r^c- Using (i) and 
(ii) one obtains 



dx2 



with boundary conditions $ ' (x)||x 



I'^l- 



0. This 



problem also can be solved in the framework of quasiclas- 



FIG. 6: Current density versus p in units of j^e — \e\k1pc/m: 
(a) T = 0, ^ = Q.^pl/m for pac = 4 and pc ~ 10, the lines 
1-3 take into account antibonding, bonding and = states, 
the line 4 includes all of them; (b) the same as in panel (a) 
for Psc ~ 2 and pc = 5; (c) T = and different chemical 
potential value p — 0.8x, 1.6x, 2.4x, 3.2x p1/2rri marked 
1-4 correspondingly (psc = 4 and pc = 10); (d) the same as 
in panel (c) for psc ~ 2 and pc = 5. 



sical approach. One has to use quasiclassical function in 
the next form (antibonding) 



(x) 



sm 



A- 



Iod^'\lk'-[T) (|x'|-xo)2 
[P-(^)'(|x|-xo)2]V4 
or (bonding) for |x| > +0 (see Appendix B and [? ]) 



(A3) 



B- 



Jodx\ k 



h 



|x'|-xo) 



[P-(^) (|x|-Xo)2]l/4 



The 



standard 



conditions 



are 



(A4) 



satisfied 



(x) = 0, (A2) if 



/, dxVA:2-(c.em/?i)2(|x'|-xo)2 



J^' dx'y^k^ ~ {ujcm/ny{\x'\ ^ xo)2 = 7r(n + 1/2 ± 1/4) 
{here n — 0,1,2..; for antibonding and "— " 

for bonding states; x^.h = ±[xo -I- hk/{ujcm)]} 
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and 



py < V2mE. For > \/2mE in 

terms of quasiclassical approach wc have four 

turning points: x^^ ^ = ^[^o ± h.k/ [ujcfh)], 

Xftj 2 = (xo ± hk/{u!c'm)) and the Born-Sommcrfeld 
quantization rule /'^"^ dx' 

leads to the spectral relation {p" < ^/2mE): 



,)2 = 7r(n 



xo)2 = 
1/2). This 



2m 



arcsm 



V hk 



TT 

+ 2 + 



hk Y 



«+^±^ (A5) 



The solution of this equation may be introduced through 
the Airy functions od first Ai and second kind Bi as 

= a±Ai[-{ci ± C2)^/^X - C3/{CI ± C2\f/^] 
+6±Bi[-(ci ± C2y/^X - C3/(C1 ± C2|)2/3], (B2) 



here "+" for a; > and "— " for a; < 0; the constant 
coefficients 



a±=Bi[-C3/{ci±C2f/% 

b± = -Ai[-Cs/{ci±C2f/% 



(B3) 



and E = huic{n + 1/2) for p^ > V2mE. The spectra 
that is definite by eq. (B5) well studied for 2D case, 
when k^ = 0. It gives usual Landau levels, when xo — > 
hk I {uicffi) ; levels splitting takes place with decreasing of 
xq; and fc — )• oo, when xq — )■ —hk/{ujcfn). Notice, that 
performed transition is valid only for k^ <^ kc- 



for asymmetric solutions, or 

a± = (ci ± C2)i/3i?z'[-C3/(ci ± C2f/% 

b± = -(Ci ± C2)1/3A*'[-C3/(C1 ± C2)2/3]. 



(B4) 



Appendix B: Solutions about magnetic field step 



Lets consider eq. (3) for |a;| -C 1. Performing lin- 



exp(2x)p^ (k^ + 2k^|x| 



earization 

« cix + C2\x\ + C3 [ci = {pcfCf I {pc. Kf- - 

C2 = pIu^ / \/ {pcK,)'^ - rri^, C3 = \/{Pci^Y - m?], we ob- 
tain 



0. 



(Bl) 



for symmetric solutions. The coefficients do not sat- 
isfy normalization condition which is not essential for the 
discussed here characteristics. One can notice that first 
derivative of symmetric solution equivalent to the sym- 
metric solution in the point a; = if Ci , C2 and C3 are the 
same in both cases. 

Similar description for the 2D or flat geometry case 
one can find in [? ]. To obtain it one has use the 
substitution x = x, ci = 0, C2 = —2xo{rhLOc/1i)'^ , 
C3 = k"^ +xl{fhujc/h)'^. 
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